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ABSTRACT 

In this paper, wc study the distribution of dark matter halos and voids using high 
resolution simulations in f{R) gravity models with the chameleon mechanism to screen 
the fifth force in dense environment. For dark matter halos, we show that the semi- 
analytic thin shell condition provides a good approximation to describe the mass and 
environmental dependence of the screening of the fifth force in halos. Due to stronger 
gravity, there are far more massive halos and large voids in ,f{R) models compared with 
ACDM models. The numbers of voids with an effective radius of 15/i~^Mpc are twice 
and four times as many as those in ACDM for f{R) models with |/flo| = 10"^ and 
10~^ respectively. This provides a new means to test the models using the upcoming 
observational data. We also find that halos inside voids are all unscreened in our 
simulations, which are ideal objects for the gravity test. 
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1 INTRODUCTION 

The biggest problem in cosmology is to explain the observed 
accelerated expansion of the universe. Within the framework 
of General Relativity (GR), the acceleration originates from 
dark energy (Copeland et al. 2006). Alternatively, a large- 
distance modification to GR may account for the late-time 
acceleration of the universe. 

It has been recognised that usually once we modify Ein- 
stein gravity on large scales, there can appear a new scalar 
degree of freedom in gravity which mediates a fifth force. 
Without a mechanism to suppress this additional force, 
most modified gravity models are excluded by stringent con- 
straints on deviations from GR on the solar system scale. 
One way to evade these constraints is to exploit a chameleon 
mechanism (Khoury & Weltman 2004). The new scalar de- 
gree of freedom couples to the energy density of matter. By 
tuning the coupling and potential for the scalar mode, it is 
possible to realise a situation that in dense environments 
such as the solar system, the scalar field has a large mass 
and it essentially does not mediate the fifth force. On the 
other hand, on cosmological scales, this scalar mode is light 
and modifies gravity significantly. 

In models with the chameleon mechanism, there ap- 
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pears environmental dependence on the properties of dark 
matter distributions. In dense environment such as in clus- 
ters, the chameleon works efficiently and the modification 
of gravity is suppressed. On the other hand, in underdense 
regions such as in voids, the chameleon mechanism does not 
work and gravity is significantly modified. As is shown in 
our previous paper (Zhao, Li & Koyama 2011b), this envi- 
ronmental dependence is a smoking gun for the modification 
of gravity in models with the chameleon mechanism. 

In this paper, we study the properties of dark matter 
halos and voids in models with the chameleon mechanism 
to reveal the environmental dependence of dark matter dis- 
tributions. We use f[R) gravity as an example and exploit 
recent results from high resolution simulations described in 
Zhao, Li & Koyama (2011a). In the f{R) gravity, the fifth 
force can enhance gravity by a factor of 1/3 but this en- 
hancement is suppressed by the chameleon mechanism in the 
overdense regions. Some other numerical simulations per- 
formed for models with chameleon mechanism are those in 
Oyaizu (2008); Oyaizu et al. (2008); Schmidt et al. (2008); 
Li & Zhao (2009); Zhao et al. (2010); Li & Zhao (2010) and 
Li et al. (2011). Although we will not study those simula- 
tions directly, we expect that the results found here will be 
true for them as well. 

The paper is organised as follows. In Section II, we sum- 
marise f{R) gravity models that we shall study in this paper. 
In Section III, we study the probability distribution function 
of the smoothed density field and show how the properties 
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of halos and voids in f{R) gravity are modified compared 
with the standard ACDM model. Section IV is devoted to 
the study of dark matter halos. We study how the differ- 
ence between dynamical and lensing masses, which arises 
due to the fifth force, depends on the mass and environment 
of halos. We find that the semi-analytic thin shell condition 
that determines the efhciency of the chameleon can well de- 
scribe those dependence found in simulations. In Section V, 
we study the underdense regions by identifying voids in our 
simulations. We study the properties of halos inside and near 
the voids. We show that the number density of large voids 
is significantly modified in /(i?) gravity models. 



2 F{R) GRAVITY AND SIMULATIONS 

The f[R) gravity, in which the Ricci scalar R in the Einstein- 
Hilbert action is generalised to a function of R, was designed 
to explain the observed cosmic acceleration without intro- 
ducing dark energy. In such theories, the structure formation 
is determined by the following equations, 

= }^aHp^+'^5R{U), (1) 
3 6 

= -y[5i?(/R)+8^G5pM], (2) 

where $ denotes the gravitational potential, fa = -^^^ 
is the scalaron, the extra scalar degree of freedom, SR — 
R — R,SpM ~ pM — PM, and the quantities with overbar 
take the background values. In GR, the gravitational po- 
tential is solely determined by the distribution of matter, 
say, V^"I> = AivGa^SpM- This is a linear Poisson equation, 
which is much easier to solve. In ,f{R), however, the scalar 
field complicates the Poisson equation, making the effective 
Newton's constant vary with the local density: in under- 
dense regions, the SR{fR) term in Eq. (1) vanishes thus two 
equations decouple, making the effective Newton's constant 
enhanced by 1/3. On the other hand, in the dense region, 
Sfa in Eq. (2) is negligible, so SR{fR) = —SttGSpm, which 
means that GR is locally restored. This is the chameleon 
mechanism making f{R) evade the stringent solar system 
tests, thus is important for the cosmological viability of the 
f{R) gravity. 

One could rewrite Eq. (1) as, 

= 47rGa^5/9eff, (3) 

where the effect of the scalar field is absorbed into the defi- 
nition of the effective energy density Sp^s- Then the dynam- 
ical mass Md (r) of a halo is defined as the mass contained 
within a radius r, inferred from the gravitational potential 
felt by a test particle at r. It is given by Mo = J a^SpcsdV, 
in which the integral is over the extension of the body. On 
the other hand, the lensing mass is the true mass of the halo, 
I.e., Ml = / a^SprndV. 

Comparing the lensing mass with the dynamical mass 
of the same halo can in principal be a easy way to test GR. 
This is because the lensing mass and the dynamical mass 
are identical in GR, but quite different in MG scenarios. To 
quantify the difference, we calculate the relative difference 
Am between Ml and Md for each halo. Am = Md/Ml — 1. 
Note that in f{R), Am ^ 1/3 (Zhao, Li & Koyama 2011b). 
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Figure 1. (Colour Online) The probability distribution of matter 
density contrast field 5 = p(x; R)/p — 1. The i5 field was filtered 
by a top-hat window with radius R = 2/i~^Mpc. In the plot we 
off'set 5 by 1 for the ease of using the logarithmic scale and the 
points plotted homogeneously in lg(l + (5). The black squares, red 
circles, green triangles and blue diamonds are from the ACDM 
simulation and f{R) simulations with Ifuo] = 10~^, 10~^, 10~* 
respectively. Each curve represents the averaged result over ten 
realisations and is normalised so that the integration of P{1 + S) 
is 1. 

The presence of the chameleon effect indicates that Eqs 
(1) and (2) are highly nonlinear, so that the system can- 
not be solved without using N-body simulations. In this 
work, we we shall use the high-resolution A^-body simula- 
tion catalogue (Zhao, Li & Koyama 2011a) for a f(R) grav- 
ity model, f{R) = aR/{l3R + -y) (Hu & Sawicki 2007) where 
a = — m^ci,/3 — C2,'y — —rn^,-iv? — HqQm and Ci,C2 are 
free parameters. The expansion rate of the universe in this 
f{R) model is determined by ci/c2, and the structure for- 
mation depends on |/i;o|, which is the value of \df/dR\ at 
2 = 0, and is proportional to ci/c^. We tune Ci/c2 to obtain 
the same expansion history as that in a ACDM model, and 
choose values for |/jjo| so that those models cannot be ruled 
out by current solar system tests. To satisfy these require- 
ments, we set ci/c2 = Qi^A/^M and simulate three models 
with |/iio| = 10~*, 10~^, 10~^. In this paper, we study the 
distribution of dark matter halos and voids in these simula- 
tions at 2 = 0. The method to identify halos is described in 
Li & Barrow (2011); Zhao, Li & Koyama (2011a). 



3 PROBABILITY DISTRIBUTION OF 
DENSITY FIELD 

In the standard cold dark matter paradigm, structures grow 
from the small inhomogeneities in the initial matter density 
field due to the pull of gravity. As a result, initial overdense 
(underdense) regions become more and more overdense (un- 
derdense). In f{R) gravity, gravity can be enhanced, so that 
the fifth force helps to pull more matter into overdense re- 
gions, and the underdense regions can be evacuated more 
efficiently. 
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In Fig. 1 we show the probabihty distribution of the 
matter density contrast field measured from our f{R) and 
ACDM simulations. This is calculated by filtering the den- 
sity field by top-hat windows with as an example radius 
R — 2/i^^Mpc centred at each cell of the simulation grid, 
and counting how many such windows fall into a given den- 
sity band. 

As shown in this figure, the fifth force can tremendously 
increase the chance of creating extremely low-density regions 
in the Universe. For example, only ~ 0.01% of the space in 
the ACDM paradigm has a density oil-\-5 = 0.05, while the 
f{R) models predict 5 (for |/flo| = 10"®), 15 (for |/flo| = 
10"^) and 30 (for |/flo| = 10"'') times as much. The effect 
becomes smaller for increasing 5, and for 1 -\- 5 ^ 0.2 — 0.3 
the probability becomes roughly the same for all models. For 
0.2 — 0.3 < 1 -|- 5 <~ 10 the fifth force actually decreases 
the probability and then for windows with 1-1-5 >~ 10 the 
fifth force makes it more likely to be found again by making 
matter cluster more strongly. Similar effects have been found 
for other models (Hellwing & Juszkiewicz 2009; Li 2011). 

The peaks of the density distribution shifts towards low 
values as |/flo| increases, which shows that the Universe in 
f{R) gravity may look emptier overall. Meanwhile, Fig. 1 
confirms that a f{R) universe will more likely host very big 
voids and very massive dark matter halos, both of which are 
rarer in a ACDM universe. We will come back to this point 
later. 

Fig. 1 also clearly shows the effect of the chameleon 
mechanism, which is known to work better for smaller j/floj 
and for high density fields (Zhao, Li & Koyama 2011a). For 
I /ho I ~ lO"'', the deviation from LCDM is suppressed for 
high density fields while there is still a sizable deviation 
in the probability distribution for under-density fields. This 
shows that the modification of gravity is more prominent 
for voids for small |/jjo|- This fact is important when we 
perform observational tests of modified gravity models with 
a realistic value of |/jjo| compatible with the solar system 
constraints. 



4 OVERDENSE REGIONS 

In general, dark matter halos reside in high-density regions, 
which form their local environment. It is well known that 
the fifth force in f{R) gravity sensitively depends on the 
environment. Thus from a theoretical point of view, it is 
very important to understand how the environment changes 
the properties of the fifth force. 

As mentioned above, ,f{R) gravity is a subclass of the 
chameleon scalar field theory, with /_h = exp(7y^<^) — 1, in 
which K — 8itG — Mp^ and (fi is the corresponding scalar 



field and 7 = 1^/2/3 is the constant coupling strength. The 
scalar field is governed by an effective potential (see, e.g., 
Li & Barrow (2007)) 



2k(1 + /h) 4 



(4) 



When the chameleon mechanism is at work, a spherical body 
will develop a thin shell, and the thickness of which is given 
by (Khoury & Weffman (2004); Li & Efstathiou (2011)), 

Afi _ (fout - yin 

R J^pinR'^ ' 



(5) 



in which R is the radius of the body, AR is the thickness 
of the shell, (y3out,V'in are the values of minimising Vefi 
inside and outside the body, respectively. Similarly, pin and 
Pout are the constant matter density inside and outside the 
body respectively. Note that only the matter inside the thin 
shell contributes to the fifth force exerting on a nearby test 
particle. From Eq. (5) it is evident that the fifth force could 
be suppressed if shell becomes thinner, which can be realised 
by the following two ways: 

(i) Increasing R and/or pin, thereby making the body (in 
the case here the dark matter halo) more massive. 

(ii) Decreasing (/Sout, which involves increasing pout or 
equivalently making the environment denser'^. 

As a result, massive halos in dense environments are strongly 
screened from the fifth force, while small halos in low-density 
environments are less screened and may experience the full 
fifth force. 

In the f{R) gravity theory, the thin-shell expression 
Eq. (5) can be translated into the following equation. 



(6) 



by using the relationship between y^(^ and fn, and the fact 
that v^iyS ~ ^ 1. The ratio between the magnitudes of 
the fifth force and gravity can be approximately estimated 
as (Li & Efstathiou 2011) 



A^ 



— X mm 



_xmm|^,l| 



(7) 



Am has a maximum value of 1 /3 as expected, and it can be 
analytically estimated by calculating AR/ R from Eq. (6) as 
follows: 

(i) Given a halo's mass and virial radius we can compute 
the average pin and therefore /ii,in; 

(ii) Pout, the environmental density, can be estimated by 
computing the average density of a sphere with a radius JJonv 
centring on the concerned halo (Li & Efstathiou 2011). Then 
/ii.out follows straightforwardly. 

For the f[R) model considered here, we have 



(l + 4|b)' 

Ha \2 



/i?.0, 



jR.,out — — -jjRO, 



(pout+4^) 



(8) 



(9) 



at z = where pout(in) = Pout{in)/p- So we get Aj\/ from 
Eq. (7), with 



AR 
R 



[1 + 4 



^ Note that the first term on the right-hand side of Eq. (4) is a 
runaway potential of ip, while the second term increases exponen- 
tially in If. So Vcs{ip) has a global minimum, which shifts towards 
smaller values of tp when 
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.(pout +4^)^ (Pin+4i^) = 



(10) 



This could be measured for each halo from the A'^-body sim- 
ulations, as is shown in Fig. 2. Note that the thin-shell con- 
dition for dark matter halos have been studied using A'^-body 
simulations by (Schmidt 2010), but here for the first time we 
have checked the environmental effects and compared with 
analytical results. 

The left panels of Fig. 2 show the dependence of Am on 
the halo mass (illustrated by the size and colour of the sym- 
bols) and the environment matter density (horizontal axis). 
The results are measured from our A'^-body simulations us- 
ing two values of fJonv, respectively 8/i~^Mpc (upper left) 
and 5/i^^Mpc (lower left). There are several interesting fea- 
tures: 

(i) Massive halos mostly reside in overdense regions, as 
could be seen from the correlation between the size of the 
symbols and the horizontal axis. This is as expected, because 
only in those regions there are enough particles to form large 
halos; 

(ii) With the same environmental density, small halos are 
less screened (have bigger Aj\/), which agrees with the anal- 
ysis above; 

(iii) For halos with comparable mass, those in overdense 
regions are more strongly screened, because of the environ- 
mental effect; 

(iv) Only very few halos reside in underdense regions, and 
those are mostly small halos. This is because most particles 
in those regions have been pulled away. Note that the halos 
in those regions are essentially unscreened; 

(v) The dependence of the results on -Rcnv is fairly weak, 
indicating that the exact value of i?cnv in the definition of 
the environment is not crucial. 

In the right panels of Fig. 2 we have shown the analyt- 
ical approximation obtained by using Eq. (10). We can find 
that the the analytical result agrees with the A^-body simu- 
lation result quite well, and once again there is no sensitive 
dependence on 7?cnv . This means that the analytical approx- 
imation Eq. (10) can well describe the nonlinear behaviour 
of the fifth force in f{R) gravity. 

In Fig. 3 we show the same results for a different 
/(i?) model, with |/flo| = 10"^ A gain we could see that 
the analytical and numerical results agree well. Note that 
here all halos but the very massive ones are unscreened. 
Figs. 2, 3 lend supports to the simple excursion-set model 
of Li & Efstathiou (2011) for studying structure formation 
in the chameleon- type scalar field models. 

Of course, the agreement between the analytical and 
simulation results is not perfect, although it is fairly good 
statistically. For |/ko| ~ 10~^, the analytic thin shell con- 
dition underestimates the fifth force especially with i?cnv = 
5/i~^Mpc. This is probably because the Chameleon mecha- 
nism becomes weaker for larger |/ho| and the environmental 
effect becomes more non-local. Thus we need a larger -Renv 
to capture the environmental effect correctly. 

Another possible reason for the slight discrepancy be- 
tween analytical and numerical results is that the dark mat- 
ter halos are not rigorously spherical, which we have as- 
sumed when applying the thin-shell condition; we have tried 
to search for the possible correlation between the ellipticity 
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Figure 5. Screening of dark matter halos (the open circles) inside 
voids in one of our |/ijol = 10~^ simulations. Only halos more 
massive than 2 -IQl^ Mq are shown for resolution considerations. 
The solid line is Aj\/ = 1/3 and we can see that these halos are 
completely unscreened. 



of the halos and Am, but we didn't find any evidence, which 
means that using the spherical thin shell condition is indeed 
a good approximation. 



5 VOIDS 

Next let us turn to the underdense regions, or voids, in the 
f{R) gravity. 

The voids are identified using VAMSUR (Voids As Merged 
Spherical Underdense Regions) code developed by Li (2011). 
The basic idea is to first find the low-density spherical re- 
gions (protovoids) and then merge them to form irregularly- 
shaped voids using a given algorithm. In this work we have 
chosen S < —0.8 as the definition of voids. 

In Fig. 4 we show the protovoids (left panel) and dark 
matter halos (right panel) identified in one of our simulation 
boxes for the model with |//?.o| = 10~^. We can see that: 

(i) Most dark matter halos (in particular the more mas- 
sive ones) distribute in regions where few protovoids can be 
identified, and vice versa, which is a trivial test of the code 
and numerical results. 

(ii) Near the voids the dark matter halos are less screened 
(denoted in blue) while far from the voids they can be well 
screened (in red). 

(iii) Halos inside voids are all unscreened as is also shown 
in Fig. 5. 

These observations agree with our expectation very well. 

As mentioned above, in the f{R) gravity the fifth force 
helps evacuate the low density regions, which results in more 
large voids than in ACDM. To see this more clearly, one 
could plot the number density of voids as a function of their 
effective volumes, and this is shown in Fig. 6 (see figure cap- 
tion for details). Here we can see the clear trend: increasing 
I /ho I, which makes the fifth force less suppressed from ear- 
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Figure 2. (Colour Online) Screening of dark matter halos as a function of the environment and halo mass for the model with |/ijo| = 10 ^. 
The horizontal axis of is ^env ) the matter overdensity of a sphere (the environment) centred at each halo, with comoving radius Rcnv as 
indicated in the panel. The vertical axis is Aj\/, the ratio between the magnitudes of the fifth force and gravity at the surface of a halo. 
Each circle represents a halo, and the halo's mass is illustrated by both the size (increasing size for increasing mass) and colour (from red 
to violet for increasing mass) of the circle. The left (right) panels are numerical (analytical) results (see text for a detailed description). 



numerical result, R = 8 Mpc/h 
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Figure 3. (Colour Online) The same as Fig. 2, but for the model with l/flo 
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Figure 4. (Colour Online) A visualisation of the distributions of voids and dark matter halos in one of our |/flo| = 10"*^ simulations. 
Each bubble represents a protovoid (left panel, where the bubble size characterises the size of the protovoid) and halo (right panel, where 
the bubble size characterises the halo's mass). The colour of a halo represents the screening of that halo, from strongly screened (red) to 
unscreened (blue). 
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Figure 6. (Colour Online) The void number density as a func- 
tion of a volume. The black squares, red circles, green triangles 
and blue diamonds are from the ACDM simulation and /(i?) sim- 
ulations with l/ijol = 10~^, 10~^, 10~^ respectively. Each curve 
is the averaged result of ten realisations. The magenta pentagons 
are results for ACDM from Colberg et al. (2005) for consistency 
check. All results arc at a = 1. 



lier times, produces more large voids. For example, a f[R) 
universe with |/ho| — 10^^ (10^*) has twice (four times) 
as many voids with effective radius 15/i~^Mpc as a ACDM 
universe does, and the relative difference in the abundance 
for larger voids is even larger. 

Very large voids and very massive dark matter halos 
are rare objects in a ACDM Universe, and we have seen that 
both of them are more abundant in the f{K) universes. This 
is the reason why cluster abundance gives the strongest con- 
straints on |/i?o| with current observations (Schmidt et al. 
2009; Ferraro et al. 2011; Lombriser et al. 2011). Compared 
with the halos, the modified gravity effect is more pro- 
nounced on the voids. This is because in f{R), gravity is 
maximumly enhanced in voids due to the presence of the 
fifth force. At 2 = 0, the model with |/i?o| = 10~^ predicts 
~30% more halos with mass ~ 5 x 10"Mq than ACDM 
(Zhao, Li & Koyama 2011a), while it predicts twice as many 
voids of size ~ 15000/t"^Mpc^. 

Due to the limitation of our simulation box size, we do 
not have voids with radius larger than ~ 15/i~^Mpc. How- 
ever, there is an abundance of such large voids observed. For 
example, using the SDSS DR7, Pan et al. (2011) identified 
about 1000 voids in the northern galactic hemisphere with 
radii > 10/i~^Mpc; the largest and median radii in their void 
catalogue are 30/i~'^Mpc and 17/i^^Mpc, respectively. Those 
voids have an edge density contrast of S < —0.85. They find 
that their observations agree quite well with the ACDM sim- 
ulations, which means that their data could place strong 
constraints on the f{R) gravity, making voids a promising 
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tool to study the physics of the accelerated cosmic expansion 
and large-scale structure formation. Of course, their voids 
are identified by looking at galaxies in the survey, and to 
make direct comparison with their data we have to generate 
galaxy catalogues in the ,f{R) gravity. We will leave this to 
future work. 



6 SUMMARY AND CONCLUSIONS 

In this paper, we studied the over- and underdense regions, 
namely, the distribution of halos and voids, in f[R) grav- 
ity simulations. By comparing the probability distribution 
function of the density contrast 5 for f{R) models with that 
for GR, we find that there are far more voids in f{R) gravity 
than that in GR. For example, the numbers of voids with 
an effective radius of 15/i~^Mpc are twice and four times as 
many as those in GR for f{R) models with |/ijo| = 10"^ and 
10~* respectively. This in principle provides a new means to 
test GR observationally using the upcoming data. We also 
find that halos near the voids are less screened and experi- 
ence stronger gravity. Especially, halos inside the voids are 
all unscreened in our simulations. This confirms the expec- 
tations that small galaxies inside voids provide us the best 
place for testing modification of gravity. 

On the other hand, the overdense regions, ie, the dis- 
tribution of dark matter halos, can provide important infor- 
mation for the GR test as well. In this work, we utilised the 
thin-shell condition developed in (Khoury & Weltman 2004) 
and (Li & Efstathiou 2011), and analytically predicted the 
fractional difference between the lensing mass and the dy- 
namical mass of dark matter halos, Aa/, as a function of 
the environment. As we found, the analytic result agrees 
very well with the simulation result, which means that the 
thin-shell condition is a good approximation for f{R) grav- 
ity. This has important applications for the semi-analytic 
halo model building for f{R) gravity, which is crucial for 
realistic constraints of f{R) models using observations. 
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